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Experiment E910 has measured slow protons and deuterons from collisions of 18 GeV/c protons 
with Be, Cu, and Au targets at the BNL AGS. These correspond to the "grey tracks" first observed 
in emulsion experiments. We report on their momentum and angular distributions and investigate 
their use in measuring the centrality of a collision, as defined by the mean number of projectile- 
nucleon interactions. The relation between the measured Ngrey and the mean number of interactions, 
V[Ngrey), is Studied using several simple models, one newly proposed, as well as the RQMD event 
generator. RQMD is shown to reproduce the Ngrey distribution, and exhibits a dependence of Ngrey 
on centrality that is similar to the behavior of the simple models. We find a strong linear dependence 
of Ngrey on u , with a constant of proportionality that varies with target. For the Au target, we 
report a relative systematic error for extracting ^{Ngrey) that lies between 10% and 20% over all 

^"^ grey • 



I. INTRODUCTION 

The use of high-energy coUisions of hadrons with nu- 
clear targets to study the space-time development of pro- 
duced particles was first suggested many years ago [Q-Q . 
Early experiments indicated that at sufficiently high en- 
ergies, the projectile will undergo on average a number of 
inelastic hadron-nucleon scatterings roughly equal to the 
mean interaction thickness, V — Aahp/crhA, with most 
particles forming well outside the target nucleus P-p). 
These data suggest that a single p-A collision can ef- 
fectively be modeled by a cascade of v proton-nucleon 
interactions, with ~1 fm formation times for produced 
particles. For reasons given below, any conflicts between 
such a cascade model and p-A data have yet to be demon- 
strated. 

The differences between a p-A collision and a p-nucleon 
cascade are especially important to discover and un- 
derstand in light of recent experiments with relativistic 
heavy-ion collisions at BNL and CERN. Here the com- 
plex hadronic physics processes that we wish to study 
in p-A form a significant background in the search for a 
QCD phase transition. The overwhelming complexity of 
A-A collisions makes it difficult to study these processes 
directly, whereas p-A collisions arc simpler and may pro- 
vide more insight. 

Many previous p-A experiments were limited by their 
inability to trigger on central collisions, those with small 
impact parameter and in which v attains the highest val- 
ues — essential to studying the effects of multiple inter- 
actions. Other experiments which could trigger on cen- 
trality were limited by low rates (low statistics) and/or 
insufficient phase space coverage for identified particles. 
However, they were able to establish a relationship be- 
tween v and a measurable observable, the number of slow 
singly charged fragments (grey tracks M ) emitted in the 
collisions. This relationship is expressed as a conditional 
probability for detecting Ngrey grey tracks given a col- 
lision in which there were v interactions, P^Ngreyli')- 
Given a distribution 7r(j/) for the number of interactions, 
the relevant quantity for measuring centrality in p-A col- 
lisions is. 



v{Ngrey) = ^ vP{N grey\v)-K{v) 



(1) 



Several forms 

have been proposed for P{Ngrey\v) [p^|-p^, yet there 
have been few systematic studies to test the validity of 
the models' assumptions and assess the accuracy of the 
extracted values oiv{Ngrey)- 

We will focus on the two models which have been 
most commonly applied to data: the Geometric Cascade 
Model (GCM) of Andersson et al. |l|], and the intra- 
nuclear cascade calculation of Hegab and Hiifner |l3^,E2[ . 
We will then present a new model which draws on el- 
ements of both. The GCM uses a normalized geomet- 



this distribution applies equally and independently to 
the distribution of grey tracks produced by each primary 
proton-nucleon scattering. This yields an analytic form 
for joint probability distribution, P{Ngrey, J^), which has, 

Ngrey OC V. (2) 

The calculation of Hegab and Hiifner performs a sum over 
the collisions of the beam and all primary struck nucle- 
ons, assuming a straight line path through the nucleus 
for all products which follows the initial impact parame- 
ter of the projectile. The mean value of this distribution 
is given approximately by, 

(3) 
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Nnrey OC V 



Despite this fundamental difference, both models and 
variations of them have successfully reproduced the 
P(Narey) distributions for a number of experiments. 
See ||,|0|jl6|-|l|l, for the GCM, and |l|,|2^ for the cas- 
cade of Hegab and Hiifner. Values of v{Ngrey) have 
been extracted for many types of experiments: emul- 
sions!^], counters [19p^, and streamer/bubble cham- 
bers |l5|l|,|0|,|| . The GCM model has also been applied 
to the Ngrey distribution from z/-Ne interactions |jl^. In 
each case, the agreement between model and data is quite 
reasonable given the simplistic nature of the models, but 
the accuracy of the v{Ngrey) extraction is undetermined. 
If the systematic errors are small compared to the range 
of V for a given target, then the analytic approaches to 
determine V(Ngrey) are justified. 

Here we present a high statistics analysis of low mo- 
menta protons from collisions of 18 GeV/c protons inci- 
dent on three nuclear targets: Be, Cu, and Au. The data 
were taken by BNL E910, a large acceptance TPC spec- 
trometer experiment with additional particle identifica- 
tion from time-of-flight (TOF) and Cerenkov (CKOV) 
detectors. To extract V{Ngrey) and assess its accuracy 
we apply several models to these data and to the dis- 
tributions produced by RQMD pSf, a cascade model for 
p-A and A-A collisions. We estimate the systematic er- 
rors inherent in the models and in the assumptions of the 
definition of Nnrm ■ 



The E910 experiment is described in Sec. ||. In Sec. IH 
we present the reduction of the data, including all cuts 
and corrections. Final results are shown in Sec. IV. 
Sec. M contains the comparisons to RQMD, and we de- 
termine the systematic errors in Sec. VI. In Sec. VH we 



ric distribution for P{Ng 



1) and assumes that 



present our conclusions. In all included figures, we will 
continue to use the term 'Wgre^" to refer to the number 
of singly charged slow fragments measured by our TPC 
in a collision, to be consistent with most of the litera- 
ture. Other commonly used terms for the grey tracks are 
"prompt protons" and "slow particles" . 



II. EXPERIMENTAL OVERVIEW 

The experimental layout for E910 is shown in Fig. 1^. 
The following discussion assumes a coordinate system 



that is right-handed Cartesian, with the beam direction 
nominally along the z-axis and the y-axis along the ver- 
tical. The time projection chamber (EOS TPC Q) has 
dimensions 96 x 75 x 154 cm, and is read out through a 
120 X 128 cathode-pad array. The TPC was placed in the 
center of the MPS magnet, which had a nominal central 
field of 0.5 T. It ran with PIO gas at atmospheric pres- 
sure with a vertical electric field of 120 V/cm. Additional 
charged particle tracking immediately downstream was 
provided by three drift chambers (DCl-3), placed near 
the end of the magnet. The drift chambers each had an 
active area of 172 x 100 cm, with 7 planes each, consisting 
of three views in x (one staggered), two in y (staggered), 
and two more views offset by ±60° from the vertical. 
The Cerenkov counter, with 139.7 x 190.5 cm aperture, 
was filled with Freon 114 and placed 4.8 m downstream 
of the target. Two mirror planes, above and below the 
vertical-mid-plane, with 48 mirrors each focused the light 
onto an equal number of phototubes at the top and bot- 
tom of the counter. The TOF wall consists of 32 coun- 
ters, each 15.2 x 178 x 4.8 cm arranged in a flat panel, 
610 X 370 X 86 cm, placed 8 m downstream and normal to 
the z-axis. Two more drift chambers (DC4-5) sat down- 
stream of the TOF wall, 9.6 and 10.1 m from the target. 
For these data, a buUseye scintillator detector was placed 
between the Cerenkov and TOF, 6.8 m from the target. 
It consisted of two scintillators, 14.6 x 30.5 cm adjacent in 
X, and behind them two more of dimensions 40.6 x 7.6 cm 
adjacent in y. 

Protons with nominal beam momenta of 6, 12, and 
18 GeV/c were normally incident on targets of Be, Cu, 
Au, and U. Only the 18 GeV/c beam and Be, Cu, and 
Au targets are included in this analysis. The targets, 
4% Be, 3% Cu and 2% Au targets were 3.9, 4.2, and 
3.4 gm/cm^ thick respectively, and were located in the 
TPC re-entrant window, 10 cm before the TPC active 
volume. Beam definition was provided by the SI and ST 
scintillators. SI was placed 3.8 m upstream of the target. 
It had dimensions 5 x 5 x 0.5 cm and was read out by 
two phototubes on opposite sides. ST, placed in front of 
the target, provided the coincidence for the trigger. It 
had dimensions 10 x 10 x 0.1 cm and was readout by a 
single phototube. Two veto counters, VI and V2 were 
used to tune the beam and to reject halo and upstream 
interactions. VI provided a 2 cm diameter circular aper- 
ture 9 cm downstream of SI, and V2 provided a 2 x 1 cm 
rounded aperture, 47 cm upstream of the target. Beam 
vectoring was achieved with two multi-wire chambers, 
A5 and A6, each with two horizontal and two vertical 
views. A5 was 10.36 m and A6 was 4.34 m upstream 
of the target. Four similar chambers with only X views, 
A1-A4, surrounded a series of six dipole magnets further 
upstream to measure the average beam momentum. Just 
upstream of A5, three beam Cerenkov counters, C1-C3, 
were placed in the beamline to reject pions and kaons 
(C2 only) in the beam. 

E910 ran in the Al secondary beam line of the AGS, 
with a typical intensity of 3 • lO^s^^. For these data 



the beam momentum was determined by A1-A4 recon- 
struction to be 17.5±0.2(sys) GeV/c. The LVLO trigger 
required a coincidence of ST and SI (which provided the 
start time for the experiment), in anti-coincidence with 
the veto and beam Cerenkov counters: 



LVLO = Si A 5r A Vi A 1/3 A Ci A C2 A C3. 



(4) 



The beam trigger definition furthermore required the ab- 
sence of a signal in Si during the preceding 1 fis. Beam 
triggers with no corresponding hit in the buUseye scintil- 
lator satisfied the interaction trigger. Final event statis- 
tics (after cuts) are given in Table H. A sample of target- 
out events were also taken. 



III. DATA REDUCTION 

All particle tracking and identification to be presented 
here comes from the TPC analysis described below. 
Time and pulse-height distributions are grouped into x- 
y clusters for each pad row using the center for the z- 
coordinate. A road-finding procedure extends the clus- 
ters along either direction to form a track. The initial 
momenta are determined from a fit to a helix, assum- 
ing a constant dipole for the field within the TPC and 
extending forward to the target. Tracks which originate 
from the target location (x^ cuts are employed) are used 
to determine the vertex. Those tracks used in the vertex 
determination are refit with fixed vertex to determine fi- 
nal momenta. All tracks must pass appropriate x^ cuts, 
have hits along ten or more pad-rows in z, and originate 
from the event vertex to be included in the Ngrey distri- 
bution. A GEANT simulation of the TPC shows the mo- 
mentum resolution for the Ngrey tracks to be dominated 
by multiple scattering, with a resolution of 15 MeV/c for 
1 GeV/c protons. A typical event is shown in Fig. El 

The TPC has good acceptance for the region forward 
of cos{9) = 0.4 and above a momentum of 100 MeV/c. 
The geometric acceptance, shown in Fig. 0, was calcu- 
lated with single track events thrown in a GEANT sim- 
ulation with multiple coulomb scattering enabled. The 
full acceptance which accounts for mis-reconstructed mo- 
mentum extends the acceptance correction to the lowest 
momentum bin, but in bins with finite geometric accep- 
tance the differences are less than 5%. 

Particles are identified in the TPC through their ion- 
ization energy loss, dE/dx, calculated using a 70% trun- 
cated mean. The distribution of dE/dx vs. momentum 
is shown is Fig. Q. The dE/dx distributions have been fit 
to the Bethe-Bloch formula with momentum dependent 
gaussian widths. This analysis docs not correct for sat- 
uration or non-linearities in the pulse-heights. Particles 
with dE/dx within 2.25i7 of that for a proton and further 
than 1.5(7 from the pion dE/dx are identified as protons. 
We require that deuterons lie within 2.25cr of dE/dx for 
a deuteron and further than 2.25cr from the proton and 
pion bands. Protons are identified up to a momentum 



of 1.2 GeV/c and deuterons up to 2.4 GeV/c. Two ad- 
ditional cuts are required to limit positron contamina- 
tion coming from photon conversions in the target (see 
Fig. p. Positive tracks within the positron dE/dx band 
are matched to the negative track with a common vertex 
which yields the smallest relative transverse momentum, 
qr = 2\pixp2\ / \ pi +P2 |. For qr < 0.037 GeV/c, 
the positive track is removed from the analysis. From an 
application of this cut to a lower momentum region we de- 
termined it to be ~50% effective for all targets. Since the 
positrons coming from tt'^s (the dominant source of pho- 
tons at these momenta) are more forward peaked than 
low momentum protons and deuterons, we furthermore 
reject positive tracks with dE/dx consistent with that 
of a positron that are forward of cos(6'/). The value 
of cos(6'/), given in Table |, was chosen separately for 
protons and deuterons for each target to minimize the 
contamination while preserving statistics. From the an- 
gular distributions of the paired-positrons we estimate 
final positron contamination to be less than 5% of the 
overall Ngrey sample for all targets. 

The buUseye interaction trigger accepts many elastic 
events and also beam events in which the beam mul- 
tiple coulomb scatters in the target and TPC. To re- 
move these we require that an event contain two or 
more charged particles emanating from the event vertex 
or a single charged particle with transverse momentum 
greater than 0.06 GeV/c and longitudinal momentum less 
than 12 GeV/c. The reconstructed vertex must lie within 
the projected x-y boundary of V2, and have a z-position 
within 2.6 cm of the centroid for Au and Cu, and 1.75 cm 
for Be. We also require at least one hit in each view of A5 
and A6 to reconstruct the beam vector for each event. All 
momenta are translated to the coordinate system aligned 
with the beam. Final Ngrey statistics are given in Table |. 

A typical energy range used to select the Ngrey tracks is 
30 < K.E. < 400 MeV |l|,|l| (0.24 <p< 0.87 GeV/c). 
The purpose of the lower bound is to reject fragmenta- 
tion products. That of the upper bound is to reduce the 
contribution from primary struck recoil protons. We ex- 
amine these cuts in light of recent multi-fragmentation 
data. The EOS collaboration has measured the proton 
fragmentation spectra in nucleus- nucleus collisions for 
1.2 GeV-A Au-I-C ||2^. The proton kinetic energy spec- 
tra show clear evidence of a kink at 30 MeV, and were 
well fit over the range 0-100 MeV by a two-component 
Maxwell-Boltzmann distribution with slope parameters 
of ~8 and ~50 MeV for the lowest multiplicity events, 
which are most similar to p-A collisions. The higher 
slope parameter is consistent with fits to spectra from 
4 GeV/c p+Pb in the range 40-150 MeV by a group 
at KEK ^^. For collisions that more closely resemble 
the data presented here, fragmentation spectra for 1- 
19 GeV/c and 80-350 GeV/c p-l-Xe have been measured, 
but only for fragments with Z > 3 p^,p8[. The fitted 
spectra in the range 10-100 MeV are consistent with the 
assertion that fragmentation spectra should appear ther- 
mal, with a temperature set by the mean Fermi momen- 



tum of the emitted fragments: T= ^^f- [g9|- Therefore, 
an appropriate lower limit for Ngrey lies near the Fernfi 
momentum, in agreement with the typical lower momen- 
tum limits for Ngrey found in the literature. 

Acceptance corrected momentum and angular distri- 
butions for protons are given in Fig. pi Distributions 
are shown only for p > 0.1 GeV/c and cos(6') > 0.3, 
where the acceptance is greater than 10%. The angu- 
lar distributions for all targets are nearly isotropic in the 
lowest bin, becoming progressively more forward peaked 
at higher momenta. The momentum distributions peak 
near 0.5 GeV/c for Au and at higher momentum for the 
lighter targets. The projections in momentum and angle 
for both protons and deuterons are shown in Fig. |6| and 

Fig.g. 

Based on these distributions, and the previous work 
in multi-fragmentation, we use a range of 0.25 < p < 
1.2 GeV/c for protons, and 0.5 < p < 2.4 GeV/c for 
deuterons for our definition of Ngrey The upper bounds 
reflect the limits of particle identification, 1.2 GeV/c for 
protons and 2.4 GeV/c for deuterons. The upper limits 
are higher than for most experiments, while the lower 
limits are comparable. We will explore the sensitivity of 
our analysis to this choice of cuts in a study of the sys- 
tematic errors presented in Sec. VI. With this definition 



of Ngrey , Fig. @ shows the corrected momentum and an- 
gular distributions for different values of Ngrey for the Au 
target. The distributions do not shift backwards of the 
TPC acceptance for large Ngrey, an effect which would 
bias our determination of v. 

The distributions are corrected for target out contribu- 
tion by subtracting the beam normalized Ngrey distribu- 
tions taken from runs with an empty target holder. After 
application of the vertex cut, this correction amounts to 
4% (12% of the Ngrey = bin) for Au, and 2% for Be 
and Cu. 

Finally, we correct for the contribution from secondary 
interactions in the target (interactions of the projectile 
with a second nucleus). The correction is performed it- 
eratively, according to Eq. g, 

Xq -(^)- 



Pn+l{Ngrey) — 6 

1 X 

2x0 



-* n \-^^ grey ) 



N„ 



Y^ Pn{i)Pn{Ngrey-t), (5) 

4=0 



where Xq is the p-A interaction length and x is the in- 
teraction thickness of the target. Convergence is rapid 
and only a few iterations are required. Corrections for 
tertiary interactions have been calculated and found to 
be negligible. The final distributions of slow protons and 
slow deuterons for all three targets are shown in Fig. g. 



IV. RESULTS 

We begin with the GCM |1G], which assumes a nor- 
malized geometric distribution of grey tracks for a single 



proton-nucleon interaction: 

P(7Vg.,,|^-l) = (l-X)X^--, X^-^ (6) 

where fi is the average measured Ngrey when v = \. Con- 
voluting V independent interactions, 



P{N,rey\l^) = ( ^'"^'l ^ j' ^ ) (1 ^ XT X""' 



X^ <"■-«. (7) 



The resulting distribution is recognizable as a negative 
binomial, where v is the standard k-parameter, and the 
mean, Ngrey (v) is given by vfj,. Taking the weighted sum 
over ly. 



^grey ^ ^T^i'^)Xgrey{l^) = 



V ■ II. 



(8) 



maximum on the number of protons that can be emitted 
from a single nucleus. The mean and dispersion for v are 
given by the probability distribution in Eq. 0, displayed 
as the open circles in Fig. |lj. 

The intra-nuclear cascade of |l3,0,p2| takes a very 



different approach in relating Ngrey to v. It assumes: 

1. all primary struck nucleons follow the initial pro- 
jectile trajectory, 

2. only secondary nucleons and a fraction (approxi- 
mately one half) of the primary protons contribute 

The full cascade calculation is solved numerically p3[ , 
but it has the feature that 77 is very nearly proportional 

to J Ngrey From this, the authors make the following 
ansatz |Il2|, 



Thus, Eq. H is satisfied, a direct consequence of the sum 
over V independent distributions. The full distribution is 
given by. 



P{Xgrey) =Y.P{NgreyW)^{v) 



(9) 



Two calculations for 7r(i/), Glauber and Hijing, are 
shown in Fig. [lC)| a for the three targets. Both calculate 
V from an optical model pO| using a value of 30 mb 
for the p-N cross-section and a Wood-Saxon distribu- 
tion of the nucleus. The Glauber calculation performs a 
numerical integration over impact parameter (b) assum- 
ing a binomial probability distribution for i^(5), where 
the mean and maximum values are given by the nuclear 
thickness. The results labeled Hijing pl| come from the 
HIJING Monte-Carlo event generator which in this con- 
text is equivalent to the LUND geometry code. The two 
distributions are similar. We use the Hijing distribution 
for all further analysis unless explicitly stated otherwise. 
Fig. 10b overlays the 7r(z/) distributions with the mea- 



sured Ngrey distributions. The similarity between them 
is what prompted the authors of ||l3] to suggest that the 
u and Ngrey distributions are correlated. 

The parameter X in Eq. is related to the mean value 
of Ngrey for a single proton-nucleon interaction, prompt- 
ing many authors to attempt to isolate the class of i' = 1 
events through multiplicity and leading particle cuts to 
determine X. In the context of the GCM, Ngrey{v — 1) 
is equal to the ratio Ngrey /v. We follow the method 
of [E9[ and allow X to be a free parameter in the fit of 
the Ngrey distributions. The results of the fits are given 
in Table O, along with the mean values, and ratio of 
Ngrey and V. Although for Au, the fitted X differs by 
2(7, the fitted values of X for the other targets are iden- 
tical to the definition of X in Eq. ||. The GCM fits are 
displayed as the dashed curves in Fig. ^. The model 
tends to fall below the data for low Ngrey , and above the 
data for high values. This is reflected in the large values 
of x^/dof. Note that the GCM distribution imposes no 



^\-^^grey) — ^ \J ^^ grey / ^^ grey • 



(10) 



Applying Eq. |I^ leads to the solid curves in Fig. |2[ which 
differ significantly from the predictions of the GCM. Fur- 
thermore, the quadratic dependence of Ngrey on v is very 
different from the linear relationship of the GCM. 

The contradictory nature of these two models led us 
to introduce another model, which allows for both a lin- 
ear and quadratic dependence of Ngrey on v, with the 
relative strengths determined by a fit to the data. The 
principal assumption is that for a given target, there ex- 
ists a relation between the mean number of grey tracks 
detected and the number of primary interactions which 
takes the form of a second degree polynomial. 



Ngrey{v) = Cq + CiV + C2V 



(11) 



We furthermore assume that the distribution is governed 
by binomial statistics; a total of Z target protons ex- 
ist which can be emitted and detected with probability 



N„. 
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(12) 



The full distribution of P{Ngrey) is again given by a 
weighted sum over n{iy) of Eq. ]|, and the coefficients 
of Eq. |ll| are derived from a fit to the data. The fit- 
ted function for this polynomial model is shown as the 
solid curve in Fig. [ll| and the coefficients are given in 
Table [II. The quadratic coefficients for both the Au 



and Cu targets were determined to be zero. For the Be 
target, the distribution does not extend far enough to al- 
low independent determination of a linear and quadratic 
coefficient. Given that in the fits to heavier targets the 
linear term is dominant and the quadratic term is negli- 
gible, we remove the quadratic component for the fits to 
the Be data. 



Fig. |Tl| shows that the polynomial model reproduces 
the data more accurately than the GCM. For a negli- 
gible quadratic term, the polynomial model differs from 
the GCM in only two respects, the presence of a constant 
term in Eq. nM and the use of binomial statistics. The 
latter is a natural choice, which conserves I^ for the nu- 
cleons, but we have not given a physical motivation for 
the constant term. To check that its inclusion does not 
alter the overall preference of the polynomial fit for a lin- 
ear dependence of Ng^ey on V, we removed the constant 
term and refit the data. The parameters are listed in Ta- 
ble [V, The resulting quadratic terms are still negligible, 
though finite, and the linear term remains the dominant 
contribution for Ngrey, even for large values of v. 

Fig. ^ gives T'{Ngrey) for all three models, and the 
dispersions for the GCM and polynomial models. The 
polynomial and GCM results are quite similar; they sel- 
dom differ by more than 15%, and never more than the 
dispersion of the GCM. In contrast, the intra-nuclear cas- 
cade differs significantly from the other two, with the 
difference increasing for the lighter targets. The joint 
distributions, P{Ngrey,i'), for p-Au are shown Fig. O. 
Here the increased dispersion for the GCM is evident, 
but otherwise the distributions again appear quite simi- 
lar. 



V. MODEL COMPARISONS 

Intra-nuclear cascade models have improved signifi- 
cantly since the work of Andersson et al. and Hegab 
and Hiifner, and are now capable of following the entire 
collision history in the context of the classical approxi- 
mations on which they are based. There are now several 
such models in the relevant energy range which have re- 
produced many features of the available data for hadron- 
hadron, hadron- nucleus, and nucleus- nucleus collisions. 
These models will ultimately provide a more accurate 
way to extract v{Ngrey)i however the large number of 
input parameters and assumptions require careful study. 
The aim of this section is to use one such model, RQMD, 
to study the implications of the GCM and polynomial 
models. This provides an additional test of the system- 
atic errors for these models. Their application to a newer 
cascade model also gives an important historical point of 
reference. 

RQMD (Relativistic Quantum-Molecular Dynamics) is 
a semi-classical cascade model for hadron-nucleus and 
nucleus-nucleus collisions Q. At AGS energies it func- 
tions primarily as a transport code for the nucleons, ex- 
cited nucleons, and produced hadrons. Particles can also 
interact through a mean field, here disabled, and string 
formation, rare at these energies. RQMD does not sim- 
ulate the nuclear fragmentation, and deuterons require 
the additional application of a coalescence calculation. 
The Ngrey count from our RQMD simulations includes 
only protons. Presumably some protons that contribute 



to Ngrey would bind with neutrons to form deuterons 
with roughly twice the momentum of the proton. These 
deuterons would then fall within the Ngrey momentum 
range for deuterons, leaving the overall Ngrey unaltered. 
A model data set of 200 K p+Au interaction events 
were generated with RQMD 2.2 running in fast cascade 
mode in the fireball approximation with all strong decays 
enforced. The RQMD output was then passed as input 
to the same GEANT simulation and track reconstruction 
used to calculate the E910 acceptance. The same momen- 
tum cuts were used to define the grey tracks, although 
proton identification was taken directly from the input. 
We did not simulate the positron contamination and no 
forward angle cuts were applied. This data set is labeled 
"RQMD E910". We also examine the full distribution 
oi Ngrey (no acceptance cuts), which includes all protons 
within the momentum range specified for Ngrey We re- 
fer to this data set as "RQMD 47r" . The Ngrey distribu- 
tion for RQMD E910 is shown in Fig. |l4[ along with the 
Ngrey distribution for the data (protons plus deuterons) . 
We see that RQMD over-predicts the middle region of 



N„. 



and under-predicts the extremes but nevertheless 



provides a reasonable description of the data. The GCM 
and polynomial fits were performed for both the RQMD 
E910 and RQMD An Ngrey distributions. The analysis 
procedure remains the same as it was for the data; the 
Hijing distribution for 7r(i/) is again used in the fit. The 
fitted functions are shown in Fig. |l5|, and the parame- 
ters are listed in Table M. The X value obtained for 
RQMD is larger than for the E910 data fits. For the 47r 
set X is larger by 35% from additional protons which 
fall outside the E910 acceptance. As with the data, the 
polynomial model gives a better description of the Ngrey 
distributions for both E910 and 47r sets. The quadratic 
coefficients are small and negative, but are not consistent 
with zero, as was seen for the data. 

The main goal in analyzing RQMD with the GCM and 
polynomial model is to compare the extracted v{Ngrey) 
values with the intrinsic v of RQMD. For RQMD, the def- 
inition of V requires some explanation. Above a certain 
energy threshold, cross-sections in RQMD are governed 
by the Additive Quark Model (AQM). A hadron which 
has one of its valence quarks assigned to a produced pion 
will have its cross-section immediately reduced by 1/3, 
to be restored after a proper time of 1 fm/c has passed. 
Therefore the distribution of the number of collisions re- 
ported for the projectile in the RQMD particle file falls 
well below distributions for Glauber and Hijing shown 
earlier. To obtain the appropriate value of u for compar- 
ison, we examine the history file and count all collisions 
suffered by particles that carry valence quarks of the pro- 
jectile. Counting for produced particles stops when the 
formation times elapse, and multiple collisions of valence 
quark-bearing particles with the same target nucleon are 
counted only once. The distribution of v calculated in 
this way is shown in Fig. Hq, along with the Glauber and 
Hijing calculations. We see that for large v, RQMD falls 
substantially below Glauber and Hijing. It is interesting 



to note that the relation of RQMD to Hijing in 7r(i^) is 
similar to its relation to the data in the Ng^ey distribution 
(see Fig. H). 

The comparison for V(Ngrey) and a[i'[Ngrey)) among 
the GCM and polynomial analyses of RQMD and their 
intrinsic values in RQMD are shown in Fig. |l^. The 
GCM and polynomial models generally differ by no more 
than one from the RQMD values in their prediction of 
T^iNgrey)- The intrinsic RQMD values are matched by 
the polynomial for the lowest Ng^ey , and by the GCM for 
Ngrey >3. The intrinsic dispersions are bounded by the 
predictions of the polynomial model below and the GCM 
above. It is also instructive to examine P{Ngrey,v) in 
slices of v, shown in Fig. |lq. The overall normalizations 
follow the behavior of Fig^e]. RQMD is above the GCM 
and polynomial distributions for small v, and below them 
for high v. The Ngrey distributions for a given v for 
RQMD are more accurately described by the polynomial 
model. 



VI. SYSTEMATIC ERRORS 

We estimate the systematic errors through a set of re- 
analyses of the polynomial model applied to p-Au data 
set with the following changes: 

historical — define Ngrey to be 0.3 < p < 1.0 GeV/c for 
protons, and 0.6 < p < 2.0 GeV/c for deuterons, 

glauber — substitute Glauber model for Hijing in the 
calculation of 7r(j/) (sec Fig. [lO| ), 



exclude — remove Nn 



bin from fit to data. 



The first two modifications are straightforward alterna- 
tives to the standard analysis. Removing the Ngrey = 
bin checks for a bias in our interaction trigger. Fig. pjjj a 
shows the rms deviations in the extracted V{Ngrey) with 
respect to the standard analysis. The historical momen- 
tum cuts show the largest discrepancy. For compari- 
son, the magnitude of the difference between intrinsic 
T'{Ngrey)^'^^^ and polynomial analysis of V{Ngrey) for 
the RQMD E910 model set is also shown in this fig- 
ure. This difference, Sv^Ngrey)^^^^^ ^ should include 
all systematic effects of this analysis in addition to sys- 
tematic errors inherent in RQMD. The dependence of 
SV{Ngrey)^'^^^^ on Ngrey should uot be taken as a true 
reflection of the behavior of the systematic errors, but 
rather an indication of their range. Note that it oscil- 
lates around the rms deviation of T'{Ngrey) for the his- 
torical analysis which is the dominant contribution to the 
systematic error. 

Fig. [l9| b shows the relative systematic error for the 
sum in quadrature of all three re-analyses. The \a sys- 
tematic error is 10-20%, peaking at Ngrey = (i^ « 3) for 
p-Au. This is significantly smaller than the dispersion, 
a{i'{Ngrey)), shown in the figure relative to V{Ngrey) for 
the standard analysis. The RQMD intrinsic difference is 



re-plotted as a relative difference to compare with our 
final estimate of the relative systematic error. 



VII. CONCLUSIONS 

We have measured the slow proton and deuteron pro- 
duction in 18 GeV/c proton collisions with three tar- 
gets. Be, Cu, and Au in the momentum range relevant 
to a determination of collision centrality. RQMD, a full 
intra-nuclear cascade model, provides reasonable agree- 
ment with the Ngrey distribution for the p-Au data. The 
simple GCM and Polynomial models are also fit to the 
data as part of the procedure to extract V{Ngrey)- The 
GCM imposes no upper bound on the number of pro- 
tons that can be emitted and therefore over-predicts the 
Ngrey distributions for all targets. Though not a per- 
fect fit (x^/dof of 1-100), the Polynomial model gives a 
better description of the data. 

We are unable to comment directly on the applicability 
of the model of Hegab and Hiifner until we can compare 
its predictions for the Ngrey distributions to data. How- 
ever, from the result of the polynomial fits we conclude 
that there is little v^ in the dependence of Ngrey on cen- 
trality, contrary to the predictions in |l^,|ljl . We cannot 
say that this contradicts results from previous experi- 
ments. The authors in |l3] compare to only one non- 
emulsion data set (where the target is known), and find 
reasonable agreement with their model. However, a later 
publication from this experiment found Ngrey to be ap- 
proximately linear in 77 |2^| , a result also obtained in [p3| . 
This evidence for a linear relation is consistent with the 
results of the polynomial fit and the central assumption 
of the GCM. The exact reason for this linear dependence 
is unknown, but we speculate that the main assumptions 
of the GCM are approximately true: an independent and 
equivalent cascade for each primary hadron-nucleon in- 
teraction. Deviations from this could be the reason for 
the presence of a finite constant term in the polynomial 
analysis. 

Our main result is the determination of centrality for a 
set of collisions from the measured Ngrey with two differ- 
ent models. The predictions of the two models differ by 
less than the predicted dispersions for most Ngrey Both 
models have been checked against a full cascade model, 
RQMD, and the intrinsic T'{Ngrey)^^^^ lies between the 
GCM and polynomial results. On the basis of the fits 
to the data, we ascribe the more accurate measure of 
T^iNgrey) to the polynomial model. Finally, we estab- 
lish a systematic error for this centrality measure that is 

10-20% of F(A^g^ey)- 
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TABLE I. Event statistics and forward angle cuts for all targets. Cuts for deuterons are in parentheses. 



Target 


Events 


cos 9f 


protons 


deuterons 


Au 


35520 


0.98 (0.97) 


56881 


10622 


Cu 


49331 


0.98 (0.96) 


45784 


6224 


Be 


100609 


0.94 (0.94) 


30622 


3366 



TABLE IL Mean values for Ngrey, i' and GCM fit parameters. 



Target 


^^ grey 


V 


X = Ngr.y/V 


ATfit 


xVdof 


Au 


1.98 


3.63 


0.353 


0.351 ± 0.001 


3.04- 10''/15 


Cu 


1.06 


2.40 


0.306 


0.306 ± 0.001 


910/12 


Be 


0.342 


1.36 


0.201 


0.201 ± 0.001 


4007/6 



TABLE in. Coefficients for polynomial fit to Ngrey- 



Target 


Co 


Cl 


C2 


x7dof 


Au 


-0.27 ± 0.02 


0.63 ± 0.01 


-0.0008 ± 0.0012 


1639/13 


Cu 


-0.17 ± 0.02 


0.51 ± 0.02 


-0.00005 ± 0.00242 


15/10 


Be 


-0.075 ± 0.008 


0.306 ± 0.006 


— 


95/5 



TABLE IV. Coefficients for polynomial fit to Ngrey with co constrained to be zero. 



Target 


Co 


Cl 


C2 


xVdof 


Au 


— 


0.439 ± 0.006 


0.019 ± 0.001 


1.06- 10''/15 


Cu 


— 


0.369 ± 0.005 


0.021 ± 0.001 


53/11 


Be 


— 


0.206 ± 0.004 


0.026 ± 0.003 


61/5 



TABLE V. Coeflicicnts for polynomial and GCM fits to Ngrey from RQMD 18 GeV/c p-Au. 



Target 


Co or X 


Cl 


C2 


x7dof 


GCM E910 


0.3605±0.0006 


— 


— 


8300/16 


GCM 4n 


0.4933±0.0006 


— 





6909/23 


Polynomial E910 


-0.136±0.009 


0.663±0.006 


-0.0089±0.0007 


80/14 


Polynomial 47r 


-0.443±0.007 


1.155±0.007 


-0.0075±0.0008 


69/21 
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FIG. 1. E910 layout 
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FIG. 2. 18 GeV/c p+Au event in EOS TPC. Reconstructed tracks are drawn for protons and deuterons which contribute to 
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FIG. 3. TPC geometric acceptance as a function of momentum and cos(^). 




momentum {GeM/c) 

FIG. 4. Ionization energy loss vs. momentum of particles with p < 2.4 GeV/c. The lines delimit the dE/dx particle 
identification cuts for protons and dcuterons described in the text. 
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FIG. 5. Acceptance corrected momentum vs. cos(9) distribution for protons. 
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FIG. 11. Log- likelihood fits to the event normalized Ngrey distributions for Be, Cu, and An targets with two models: The 
Geometric Cascade Model (dashed lines) and the Polynomial Model (solid lines). 
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FIG. 18. Event normalized Ngrey distributions for a given value of v for RQMD tracks in the E910 acceptance, overlayed 
with the predictions of the two models — dashed for the GCM, and dotted for the Polynomial. 
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where the prime indicates an alternate set of systematic 



cuts, or the intrinsic value from RQMD. b) the relative systematic error (shaded) for the combined one standard deviation 
errors from the three re-analyses described in the text. For comparison the RQMD relative intrinsic difference (triangles) and 
the polynomial model relative dispersion (unshaded) are also shown. 
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